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ABSTRACT 

Context. Magnetic cataclysmic variables are close binary systems containing a strongly magnetized white dwarf that accretes matter 
coming from an M-dwarf companion. The high magnetic field strength leads to the formation of an accretion column instead of an 
accretion disk. High-energy radiation coming from those objects is emitted from the column close to the white dwarf photosphere at 
the impact region. Its properties depend on the characteristics of the white dwarf and an accurate accretion column model allows the 
properties of the binary system to be inferred, such as the white dwarf mass, its magnetic field, and the accretion rate. 

Aims. We study the temporal and spectral behaviour of the accretion region and use the tools we developed to accurately connect the 
simulation results to the X-ray and optical astronomical observations. 

Methods. The radiation hydrodynamics code Hades was adapted to simulate this specific accretion phenomena. Classical approaches 
were used to model the radiative losses of the two main radiative processes: bremsstrahlung and cyclotron. Synthetic light curves and 
X-ray spectra were extracted from numerical simulations. A fast Fourier analysis was performed on the simulated light curves. The 
oscillation frequencies and amplitudes in the X-ray and optical domains are studied to compare those numerical results to observational 
ones. Different dimensional formulae were developed to complete the numerical evaluations. 

Results. The complete characterization of the emitting region is described for the two main radiative regimes: when only the 
bremsstrahlung losses and when both cyclotron and bremsstrahlung losses are considered. The effect of the non-linear cooling in¬ 
stability regime on the accretion column behaviour is analysed. Variation in luminosity on short timescales (~ 1 s quasi-periodic 
oscillations) is an expected consequence of this specific dynamic. The importance of secondary shock instability on the quasi-periodic 
oscillation phenomenon is discussed. The stabilization effect of the cyclotron process is confirmed by our numerical simulations, as 
well as the power distribution in the various modes of oscillation. 

Key words. Physical data and processes: accretion, instabilities, shock waves - Stars: white dwarfs, cataclysmic variables - X-rays: 
binaries 


1. Introduction 


Magnetic cataclysmic variables are close binary systems con¬ 
taining a highly magnetized (B ~ 10 - 200 MG) white dwarf 
that accretes matter from a low-mass companion. For a subclass 
of those objects, called AM Her stars or polars, the magnetic 
field of the white dwarf is strong enough to lock the system 
into a synchronous rotation and to prevent the formation of an 
accretion disk ( [Cropper |1990[[^rner [1995] ). Matter leaves the 
companion, follows a ballistic trajectory, and is then captured 
by the magnetic field lines of the white dwarf. Its falls down 
to the magnetic pole of the compact star at supersonic free-fall 
velocity (vff ~ 5 - 10 x 10^ cm/s), forming an accretion column 
structure. 


The collision of the infalling flow with the white dwarf 
photosphere leads to forming a reverse shock, which propagates 
upstream and heats the plasma to temperature Ts - 10-50 keV. 
These temperatures are due to the high accretion velocity 
(Ts oc v^), which is directly linked to the mass (and radius) of 
the white dwarf. Those extreme temperatures, as well as the 
strong magnetic field, imply that the post-shock region radiates 
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through several radiative processes from the optical to X-ray 
domains. The radiation arising from these binary systems comes 
mainly from the accretion region. As a consequence, these 
systems can be used as a perfect probe for testing accretion 
models in high-energy regimes. The energy losses shape the 
density and temperature profiles in the post-shock region and 
induce a strong gradient of density, temperature, and velocity 
near the white dwarf surface. Since the radiation properties 
depend on the mass of the white dwarf, modelling the emitting 
region can be used to determine the physical properties of th e 
compact object (see jKing & Lasota |1979[|Cropper et al. |1998| ). 
Numerous studies have shown, in fact, that modelling of the 
post-shock region physics s trongly affects any determination of 
the white dwarf properties (jlshida et al. |199f| jWu et al. 1 19951 


Yuasa et al. |2010t[Hayashi & Ishida |2014|). 


A stationary regime is frequently a ssumed to be established 
to help in interpreting observations ( [Fisher & Beuermann 
[200 1[ ). However, radiative instabilities that can potentially 
develop in the post-shock region can induce dynamical ef¬ 
fects and lead to the unsteady behaviour of the post-shock 
accretion column (hereafter PS AC). Indeed, the efficiency of 
the bremsstrahlung process can trigger the cooling instability. 
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The shock thus oscillates around the stationary position. First 
numerical studies of such accretion flows have highlighted the 
presence of the cooling instability in th e column ([Langer et al 


1981 ) supported by theoretical studies ( [Chevalier & Imamura 


1982 ). A comprehensive review of the different stationary and 
instationary models was given later by jWu | ( |2000| ). 


The oscillation of the front shock leads to a variation in 
the luminosity due to the induced variation of the emission 
volume, as well as of the conditions of density and temperatur e 
inside the column ( jlmamura et al. |[T991[ jWu et al~||1992j ). 
The characteristic timescale of this phenomenon is the cool¬ 
ing timescale associated to the bremsstrahlung emission, i.e. 
^brem ~ 1 s. This timescale value is consistent with the periods 
of the quasi-periodical oscillations (QPOs) observed in the 


et al. [[19831 [Bonnet-Bidaud et al. [[19851 [Larsson [[1987 

■fT^ 

Bonnet-Bidaud et al. [[1996[). These rapid oscillations are char- 

acterized by an amplitude of 1-5% rm 
From the model of cooling instability 
optical oscillations is theoretically ex 
been extensively searched in X-ray 1 
are presented in a related paper (Bom 
hereafter referred as Paper I). 

In this context, various analytical ai 
been performed to determine the sensi 
physical parameters, such as the reve 
(Strickland & Blondin |1995[[Ramach( 

the influence of different possible phys 
in the post-shock region, like two-temj 

s and a period of 1-2.5 s. 

, an X-ray counterpart of 
pected. Thus QPOs have 
ight curves. New results 
let-Bidaud et al. (2014J, 

ad numerical studies have 
tivity of the instability to 
Tse shock Mach number 
mdran & Smith |2006[) or 
ical processes taking part 

3erature effects ( 

Imamura[ 

[et al. [[19961 [Saxton & Wu [[19991 

Saxton et al. 

2005|i or 


transverse magnetic field ( Toth & Draine 1 1993 [ [Ramachandranj 


|& Smith ~||2005j ). A few studies have also described the effect 
of the coupling between several cooling processes ( jWu et ak] 
" T992t [Saxton et al. ||1998] ). In particular, recent astronomical 


observations (see Paper I) show that several polars are in 
regimes where either the bremsstrahlung dominates the losses or 
both bremsstrahlung and cyclotron have to be taken into account. 

The two main conclusions of those previous works are the 
following: Firstly, the cyclotron losses damp the oscillation. 
Thus, a critical magnetic field amplitude can be found for which 
the oscillation disappears after a few periods. To sustain an 
oscillation for a longer time, a noisy accretion rate has been 
invoked ( [Wolff et al. |lWl[Wood et al. |1992l[WiI et al. |1992[ ). 
Secondly, the oscillations of the emission in the X-ray domain 
due to the bremsstrahlung process and in the optical domain 
due to the cyclotron one are expected to present the same 
frequencies and the same order of magnitude of the relative 
amplitudes ( [Wu et al~~[[1992[ ). To complement those results 
and to better understand the physics of this specific regime, 
we present new numerical simulations to study the unsteady 
behaviour of the accretion column. Those simulations are 
compared to observational results (see Paper I). 

In Section the standard model is reviewed in order to 
build the parameter space and the relevant accretion regimes 
of the column as a function of the three leading parameters 
(specific accretion rate, white dwarf mass, and magnetic field). 
In Paper I, four parameters are considered because the specific 
accretion rate is not reached directly by observations but 
through the combination of the total accretion rate and the 
column cross-section. In Section the general behaviour of a 


column in an unsteady state is studied with particular attention 
to the additional instabilities that can develop in the post-shock 
region. Then in Section the dependency of the observational 
QPO parameters on the three leading parameters is developed 
further and discussed. 


2. Models of radiation and magnetic processes in 
high-energy accretion regions 

2.1. Brief review of the accretion column standard model 

In the standard model of accretion column, the matter coming 
from the companion is trapped by the gravitational field of 
the compact object and falls at free-fall velocity, as given by 
Vff = (2 GMwd/^wd)^^^, onto the magnetic poles of the white 
dwarf, with G, Mwd, and /?wd being the gravitational constant, 
the white dwarf mass, and radius, respectively. It should be 
noted that the matter is actually captured at the Alfven radius, 
Ra, so the free fall velocity should be corrected by a factor 
(1 - 7 ?wd/^a)- However, for polars, the Alfven radius is much 
larger than the white dwarf radius so this effect has been 
neglected. Considering the Nauenberg mass-radius relation 
( [Nauenber^[1972[ ), the free-fall velocity only depends on the 
white dwarf mass. Typically, this velocity is Vff ~ 8.75 x 10^ 
cm/s for a 1 Mq white dwarf. The incoming flow is thus highly 
supersonic with a Mach number of Al ~ (10 - 100). 

Owing to the magnetic collimation, matter spreads over a 
small fraction of the white dwarf surface, denoted / and defined 
by / = 5 ' /(47rR^j^) where S is the area of the column cross- 
section. This quantity is not well constrained by observations. 
As an example, thi s quantity has be en evaluated in the range 


/ ~ 10“^-3.10“^ by [Imamura [ ( [1984| ) for a sample of polars. The 
impact of the accreting flow onto the white dwarf surface leads 
to the formation of a reverse shock, which counter-propagates in 
the incoming plasma. The shock transforms the kinetic energy 
into thermal energy, and the PSAC is strongly heated to high- 
temperature Ts. Using the Rankine-Hugoniot conditions and as¬ 
suming that the accreted medium is described by the perfect gas 
approximation, the post-shock temperature scales as v 
given by 


and is 


1 ) 2 , , ,, 

Ts = -—;- —rv^ = 1.1 X 10 K 

kB (r+1)^ 


Vff 


10^ cm/s 


( 1 ) 


with yu, mn, and y being the mean molecular weight, the 
hydrogen mass, the Boltzmann constant, and the adiabatic 
index, respectively. We use p = 0.5 (fully ionized hydrogen) 
throughout the paper and y = 515. 

At such high temperatures, matter radiates strongly through 
optically thin bremsstrahlung, which dominates the hard X-ray 
emission ( [Frank et al. [[2002] ). Moreover, for the highest mag¬ 
netized polars, the energy can also be transported efficiently by 
the cyclotron process ( [Saxton et al. 11998] ), which emits mainly 
in the optical domain for the range of polar magnetic fields. As 
infalling material cools and slows down, it also gets denser with 
a sharp gradient near the surface. Thanks to the energy losses 
by radiation, the reverse shock position eventually stabilizes. 
Consequently, when using the steady state approximation and a 
cooling dominated by bremsstrahlung losses, the height of post¬ 
shock accretion region can be estimated and is given by [Wu et[ 
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|an ( [T994l l 


as 


3 X 10^ cm 


m 

-1 

Mwd 

3/2 

^WD 

_ 1 g/cm^/s _ 


0.5 M©_ 


. 10^ cm. 


-3/2 


( 2 ) 


with m = poVff the specific accretion rate. A more complex but 
exact analytic expression can be found in |Falize et al. | ( [2009] ). 
With the Nauenberg relation, the relation ^ can be written 
as a function of m and Mwd- Given the order of magnitude of 
the accretion shock height, the local white dwarf gravitational 
field g does not vary enough to strongly modify the post-shock 
structure (gVs « Vff). This approximation is valid as long as the 
white dwarf mass is below 1 M© ( [Cropper et al. |1999| ). 


On the whole, three key parameters determine the post¬ 
shock structure: the white dwarf mass, Mwd, determining the 
free-fall velocity and the post-shock temperature; the specific 
accretion rate, m, determining the post-shock density; and 
the local magnetic field strength, B, which determines the 
relative importance of cyclotron over bremsstrahlung losses. By 
comparing the characteristic timescales of the different physical 
process es at the front shock, the possible post-shock regimes are 
defined ( [Lamb & Masters |1979[|Langer et al.|1982| ). 

Four main regimes are obtained and can be explicitly defined in 
the (Mwd, B) diagram. Figure 1 displays these regimes for a 
0.8 M© white dwarf in the (m, B) plane. They are delimited by 
lines corresponding to an equal balance between two processes 
such as bremsstrahlung-cyclotron, electron-ion collisions, and 
ion-ion collisions. In the following, the characteristic timescales 
of these processes are recalled and used to derive the (m, B) 
relations corresponding to the limits of the regimes. 


Among all the radiative effects occurring in dense and hot 
magnetized plasmas, the dominant processes that shape the 
PSAC structure are the bremsstrahlung and cyclotron emissions 
( [Lamb & Masters |1979[|Wu |2000| ). Through the estimate of op- 
tical depth of the po st-shock medium to bremsstrahlung photons 
( Frank et al. |2QQ2| ), it can be deduced that radiation escapes the 
matter freely. Thus the loss of energy can be modelled by a cool¬ 
ing term that can be analytically expressed (jRybicki & Lightmanj 
[[T979l )as 

" (4^ ^ = Aop^T‘/2 , (3) 

where gs, nie, e, Z, Ue, nt, and p are the Gaunt factor, the 
electron mass, the electron charge, the mean charge num¬ 
ber, the electronic and ionic density, and the gas density, 
respectively. Assuming that the plasma consists of fully 
ionized hydrogen, the parameter Aq in 0 is estimated to 
Ao ~ 5 X 10^^ erg.cm^.g“^.K“^/^.s“^. According to its definition 


jRybicki & Lightman 

\\919 

), the Gaunt factor is taken as 

f ^ 1 (cf. [Saxton & Wu 




by 

^brem 


The bremsstrahlung cooling timescale, noted ^brem is given 
pe 


^brem 

^brem ~ 4 X 10 S 


m 


Ig/cmVs/ \ 10^ cm/s 


Vff 


(4) 


with e the internal energy per unit of mass. Parameters in rela¬ 
tion 0 have been normalized to the typical values encountered 


in polars. Thus, for a white dwarf of 0.8 M© with m ~ 1 g/cm^s, 
the cooling timescale associated to the bremsstrahlung process 

is ^brem ~ 1 S. 


When the local magnetic field close to the white dwarf sur¬ 
face is strong enough, the cyclotron process can also be an effi¬ 
cient radiative process for cooling the medium. Depending on 
the radiation frequency, the plasma is optically thick or opti¬ 
cally thin to photons which strongly complexifies the modelling 
( jLanger et al.|1982[|Kalomeni et al. [2005] ). Theoretical simpli¬ 
fications and dimensional considerations allow the cyclotronic 
transport to be reduced to a cooling function ( jLanger et al.|1982t 


Saxton |[1999| ), expressed as a function of the column cross- 


section area (S), the magnetic field, and the local thermodynam¬ 
ical conditions as 


A.,., = L2xlO«erg/cm>/s(^5^) 

,{ p r( ^ f 

\4.10-8 g/cmV [iO^kJ ■ 


(5) 


For a given accretion column surface and a given white dwarf 
magnetic field, <0 is reduced to 


A, 


cycl 




( 6 ) 


where Ao,c is thus a constant. 


As for the bremsstrahlung, the cyclotron cooling timescale is 
given by 


Lycl 

Lycl 


pe 


K 


cycl 


^ 2 X 10^ s ^ 


1 g/cmVsj 


/ 5 

' ^ y 

\ 10^^ cm^ / ' 

GO mg; 

\ 17/20/ 

\-77/20 

\ / Vff 

1 


-57/20 


^ 10^ cm/s j 


(7) 


As expected, the cyclotronic effects are more efficient when the 
magnetic field is stronger and the plasma is hotter. The optically 
thick property of the medium can be retrieved from ([7]), since 
tcyci increases with the surface and the density of the column. 
Indeed, the wider and denser the column, the longer it takes the 
radiation to escape the column. 


To determine the prevailing mechanism between 
bremsstrahlung and cyclotron processes, the two cooling 
timescales computed right below the sh ock front are com pared 
through a dimensionless ratio, noted 6s ( jWu et al~||1994| ), and 
given by 




^brem 

= 2 X lO"'^ 

/ 5 ' 

-17/40 

Lycl 

\10^^ cm^y 

) 


B / 

m '' 

-37/20 


10 MG/ \lg/cmVs/ 


Vff 


10^ cm/s 


117/20 


.( 8 ) 


The bremsstrahlung is the dominating process when 6s < 1 near 
the front shock, while the cyclotron is the dominant one when 
6s > 1. Since 6s oc the first regime is obtained 

for high accretion rates and weak magnetic fields {Regime 1 in 
FigQ, and the second regime corresponds to strong magnetic 
fields and low accretion rates {Regime 2 in Fig[^. 
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From Eq. ([5]) a critical magnetic field depending onm,S, and 
Vff can be obtained at the limit = I, which separates the two 
regimes and is computed as 




es=l 


200 MG 


m 


1 g/cm /s 
S \17/114 


37/57 / 

110^ cm/s 


Vff 


- 39/19 


(lO^^ cm^) 


(9) 


With the free-fall velocity expression and the white dwarf mass- 
radius relation, can be written as only a function of m, 
S , and Mwd • H should be emphasized that the bremsstrahlung 
losses increase toward s the surface, while the cyclotron losses 
decrease (see section |4.2| ). In fact the radiative losses lead 
to the cooling of the plasma and its densification. This evo¬ 
lution implies an increase in the cyclotron cooling timescale 
(^cyci ^ pi7/20 j-3/2^ ^ decrease in the bremsstrahlung scale 

(^brem P~^ Vt). Thus, even if the cyclotron dominates at the 
front shock, the bremsstrahlung will become the dominant pro¬ 
cess lower in the column, inside the high-energy region (see 
Sect.[42l). 



Fig. 1. Different regimes of accretion as a function of the accretion 
rate and the magnetic field strength, considering a white dwarf of 
M = 0.8 Mq and an accretion section of 5 = 10^^ cm^. 


The one-temperature approximation is valid as long as the 
Coulombian collisions between ions and electrons are sufficient 
to insure that the two populations have the same temperature de¬ 
spite some energy losses by radiation. When it is not the case, the 
electrons and ions are thermally decoupled and two-temperature 
effects appear. Th e characteristic timescale of electron-ion colli¬ 
sions is given by ( [Spitzer |1965| ) 


= 6.22 X 10 ^ s 


InAcV 

20 j 


m 


1 g/cm7s 


\ 4 

Vff I 
10^ cm/s / 


, ( 10 ) 


where In Ac is the Coulomb logarithm. These effects are signif¬ 
icant when is greater than the timescale of the most efficient 
cooling process. In the accretion column case, it happens when 
4i ~ ^cyci {Regime 3 in Fig[^. A critical magnetic field can be 
estimated, as for the previous processes, and is given by 




In Ac 


1 g/cm7s 


Vff y 

10^ cm/s / 


\ 10^^ cm^ / 


( 11 ) 
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with^o = 1.9 X 10^ MG. 


Finally for extremely low accretion rates, as observed at 
some epochs in most polars ( [Cropper |1990| ), the flow becomes 
collisionless (ta > t^ych where ta is the ion-ion collision 
timescale) and no shock is formed at the impact {regime 4 in 
Fig[TJ. This regime has been identified observationally among 
difiTerent polars referred to as FARPs (low accretion rate polars, 
se ejSchwope et al.' ( |2QQ2|)) . By using the expression of ta given 
in jKylafis & Famb pl982| ), this specific regime turns out to be 
crucial when the magnetic field is above a critical value that has 
the same dependency on the accretion rate, free-fall velocity, 
and accretion surface as in ( pTj ) but with a different value of the 
constant = 8.5 x 10^ MG. 


Considering more typical conditions commonly found in po¬ 
lars (see Paper I), we focus our numerical studies on the first two 
regimes {Regimes 1 and 2), which are the only ones considered 
in the following sections. 


2.2. Radiation hydrodynannic numerical code HADES 

The post-shock region is structured by the cooling processes, 
thus the position of the reverse shock can be evaluated using 
Equation For typical values of density and velocity of the 
incoming fio^ Vs ~ 3 x 10^ cm considering bremsstrahlung 
cooling alon^ As a consequence, the curvature of the magnetic 
field line close to the pole can be neglected, and the accretion 
column can be assimilated to a cylinder. Moreover, considering 
the spatial extension of the column, the magnetic field gradient 
can be neglected along the post-shock region, and the terms 
linked to the magnetic field in the ideal MHD description 
vanished so an hydrodynamic al description is sufficient to 
model the base of the column ([Cropper jj 19901 [Warner [[1995[ 
[Wu |2000[ ). 

The code Hades is d edicated to the resoluti on of radiative hy¬ 
drodynamics equations ( [Michaut et al. |2011| >. Considering our 
system, we treat radiative losses as a source term in the energy 
conservation equation 


+ V.(pv) = 0 , (12) 

d 

—(pv) -b V. (pv 0 V -b P./) = 0 , (13) 

ot ^ ' 

^ |pe + + V. ^pe\ + Ipvv^ + Pvj = -A(p, T), (14) 

where t, p, P, v, e, A(p, P), and 7 are the time, density, pressure, 
velocity, internal energy per unit of mass, cooling function of 
the system, and the identity matrix, respectively. Equations 
{Hi, ([I3]l, and ([T§ come from the conservation of the mass, 
impulsion, and energy of the system, respectively. 


The perfect gas approximation is assumed to close the 
system of equations. The code is designed to solve high-Mach 
number flows (Al >10) thanks to a MUSCF-Hancock scheme 
( [Van Feer |1979[) us ing a HFFE solver to compute the Riemann 
problem ( [Einfeldt [[1988[ ). The cooling function A is expressed 
as a sum of the power laws of density and of pressure, thus 
modelling the various cooling processes at stake in the shocked 

^ Adding the cyclotron contribution would induce a shorter cooling 
timescale, hence a shorter column. 
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medium ( [Busschaert et al. ||2013a| ). In this numerical model, 
the initial conditions are the mass density (po), the free-fall 
velocity (vff), and the temperature of the infalling flow, respec¬ 
tively. In the presented simulations, the infalling flow velocity 
evolves from 3 x 10^ cm/s (Mwd = 0.4 M©) to 5.5 x 10^ cm/s 
(Mwd = 0.8 Mq), and the infalling flow temperature has been 
set at r = 10^ K to ensure a Mach number in the flow higher 
than ten to obtain a strong shock but also to ensure numerical 
robustness since the code cannot sustain Mach numbers that are 
too high. 


Moreover, a speciflc boundary condition has been developed 
to mimic the accretion onto the white dwarf at the bottom of 
the column. Indeed, since the cooling processes stop when the 
plasma is sufficiently cooled, a layer of cold and dense mate¬ 
rial accumulates at the white dwarf surface. On the astrophysi- 
cal scale, this matter is assimilated by the star, so to avoid this 
unphysical accumulation, the mass flow is conserved along the 
column. Therefore at the bottom of the column, the boundary 
condition satisfies p^ax x ^min = Po x Vff. The parameter Vmin is 
determined so that the velocity contrast does not exceed three 
orders of magnitude, insuring the calculation’s robustness. We 
verified that the minimum temperature reached resembles tem¬ 
perature at which the cooling is no longer efficient. More details 
can be found in jBusschaert et al.] ( |2013a| ) concerning the imple¬ 
mentation and validation of this boundary condition. 


height variation induces a modification of the post-shock volume 
of the hot and highly radiating material, as well as a modification 
of the post-shock temperature. As a consequence, the integrated 
luminosity is also modified. 



Fig. 2. Evolution of the front shock position measured from the white 
dwarf surface (top blue line) and the X-ray (0.5 - 10 keV) luminosity 
(bottom red line), for a PS AC with an accretion rate of m = 10 g/cm^s 
onto a 0.8 M© white dwarf and a pure bremsstrahlung cooling. 


3. Non-stationary phenomena in the post-shock 
region 

The structure and dynamics of the PSAC can be strongly modi¬ 
fied by unsteady phenomena. In particular the cooling instability 
( [Chevalier & Imamura |1982{ [Mignone |2005[ ) can develop, and 
the post-shock matter can radiate so strongly that it again be¬ 
comes supersonic close to the white dwarf surface, which leads 
to the formation of a secondary shock ( [Falle |198T] ). These un¬ 
steady phenomena imply that no stationary regime can be estab¬ 
lished. These time variabilities are used to explain observations 
of QPOs in the optical luminosity curves of some polars. The 
effects of these two mechanisms onto astronomical observations 
are discussed. Our aim is to quantify the impact of these effects 
to evaluate the relevancy of the stationary model. These condi¬ 
tions are critical since they are used to determine the white dwarf 
mass ( [Wu et al. 119951 ). 


3.1. The cooling instability 

Thanks to the bremsstrahlung cooling, the post-shock region 
can d eveloped the cooling instability ( [Chevalier & Imamura [ 
[1982[ ). Indeed, according to the Field criterium ( [Field |1965[ ), a 
radiating fluid characterized by a cooling function in the form 
A cx is instable if JS < 1. The linear analysis has been 
the subject of various studies where the stationary solution is 
perturbated. The non-linear regime can be explored in details 
with numerical hydrodynamical simulations ( [Imamura et al. [ 
[1984[ [Mignone |2()05] ). 


Figure shows the dynamics of the accretion shock simu¬ 
lated using the Hades code with a cooling function A oc P^T^'\ 
relevant for a pure bremsstrahlung cooling (see Eq. 4). The po¬ 
sition of the front shock with time is plotted, showing a quasi- 
periodic oscillation with a period of the order of tbrem (see Eq. 5), 
in accordance with previous works (Mignone 2005). This shock 


The X-ray luminosity and spectrum are computed using the 
local bremsstrahlung emissivity combined with the density and 
temperature profiles extracted from the simulations. To compare 
with observations, the integrated X-ray luminosity is evaluated 
by integrating the radiative losses both over the volume of the 
column and over photon energies between 0.5 and 10 keV that 
correspond to the XMM-Newton detector range used in Paper I. 
The resulting integrated X-ray light curve is also shown in Fig|^ 
and clearly displays oscillations at the same main frequency than 
those observed for the shock height but with a more complex 
structure. 


3.2. Formation of secondary shock: Falle criterion 

The careful examination of the temperature and density profiles 
of the postshock region and their evolution with time clearly 
show the formation of a secondary shock (see Fig. [^. This 
secondary shock, which has been somewhat overlooked by 
previous works, is associated with the development of the 
catastrophic instability, or Falle instability, that can take place 
for a post-shock region cooled mainly by the bremsstrahlung 
process. This instability was first evoked by [Falle [ ( [1981[ ) in the 
context of the fllamentation of supemovae remnants. Secondary 
shocks in the post-shock region can arise since right below the 
main shock, the fluid velocity is subsonic, but owing to the 
cooling, the flow can again become supersonic close to the 
white dwarf surface. This kind of phenomenon has already been 


noted in several studies concerning accretion shocks (Blondin[ 

& Cioffi 119891 

Strickland & Blondin 1 19951 [Wolff et al. |1989l 

Mignone [[2005 

1, but its impact onto the main shock dynamics 


or onto the observables has not been discussed yet. 


To determine the favourable conditions for developing the 
catastrophic cooling instability, a local characteristic cooling 
length, Lcooh is defined as the product of the sound velocity, Cs, 
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and the cooling time scale, ^cooi as 


^cool = Cs X fcool oc V' “ > 


(15) 


where a general cooling function A = Aop^T^ has been consid¬ 
ered. Then, an acoustic wave propagates on a length of the order 
of Lcooi before its dynamics is strongly impacted by the cooling 
processes. It can be shown that to develop the catastrophic cool¬ 
ing instability, the local cooling length has to be smaller than the 
local distance from the surface of the white dwarf ( |Falle |19^ 


^cool 


< 1 , 


( 16 ) 


and considering a cooling function characterized by a parameter 
a = 2, SiS for the bremsstrahlung process, then a second criterion 
is for the second parameter to satisfy ( |Falle |19^ 

yS < 3/2 . (17) 


The evolutions of density and temperature in the post-shock 
region are presented at four representative times in Fig. Con¬ 
sidering a system cooled by the bremsstrahlung process, the con¬ 
dition ^rj\ is satisfied since in this case p = 1/2. The first con¬ 
dition is more complicated to check out. We calculated the ra¬ 
tio Lcooi!X over the column, which decreases from the shock to¬ 
wards the surface. When the secondary shock appears, the ratio 
tends towards 1 by superior value close to the surface. To ensure 
that the structure observed is truly a shock, we computed the 
Mach number associated with this structure over an oscillation 
(see Fig.[^ and checked that the value is around 1.2-2. Thus, 
assuming a system cooled only by bremsstrahlung, a secondary 
shock appears in the post-shock region (see Fig[^ t = 6.4 s). The 
secondary shock is formed near the surface of the white dwarf 
and progresses upstream up to the main reverse shock (see Fig|^ 
t = 6.43 s) until the two shocks collide. 


The positions of the main and secondary shocks have 
been tracked and are presented in Fig. together with the 
Mach number value for the two shocks. The interaction of the 
two shocks across the oscillation can be fully described. At 
maximum extension, the main shock progressively slows down 
and then falls back towards the white dwarf surface due to 
the cooling instability (Fig. a’). Meanwhile, the secondary 
shock develops from the white dwarf surface (Fig. a) and 
propagates upwards until it collides with the main shock (Fig.|^ 
b). After both shocks have collided, the main shock propagates 
upstream, and a rarefaction wave ( IZeTdovitch & Raizer ||2002| ) 
propagates through the shocked medium towards the white 
dwarf surface, as expected (see Figj^ t = 6.48 s and Fig.|4] c’). 
In the meantime, a new secondary shock has formed close to 
the white dwarf surface (Fig. c). On the temperature profile 
(see Fig. [^, both the rarefaction wave and the secondary shock 
are clearly visible dX t = 6.48 s near x ~ 5-6x10^ cm and 
X ~ 2 X 10^ cm, respectively. As the rarefaction wave encounters 
the newly formed secondary shock (see Fig. d), the shock 
accelerates further and then collides again with the main shock. 
Then this cyclic behaviour starts again. The secondary shock 
always appears close to the white dwarf surface when the 
extension of the PSAC is close to its maximum. After each 
collision between the two shocks, only the main shock remains 
and is pushed away from the white dwarf surface. 


Secondary shocks are also present in previous studies, 
which considered at least four different boundary conditions 
(see Fig. 5 in |Mignone~| ( [2005| )), though such shocks were 
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Fig. 3. Density (top) and temperature (bottom) profiles of the column 
at successive times across the oscillation (shown in the inset) for an ac¬ 
cretion rate of m = 10 g/cm^s, a white dwarf of 0.8 M©, and a pure 
bremsstrahlung cooling. Starting after t = 6.35 s, a secondary shock 
clearly forms at the basis of the accretion column and propagates up¬ 
stream while the main shock falls towards the white dwarf surface. 


not fully discussed in this work. To ensure that this physical 
phenomena truly exists, we verified that its presence and this 
cyclic behaviour are independent of boundary conditions at the 
basis of the accretion column and initial conditions. 

To study the development of the catastrophic cooling 
instabilities, several cooling functions have been tested with the 
parameters a = 2 and various p values. For values of p higher 
than 0.7, this instability does not develop. Indeed, even if the 
second criterion (see Eq.(p^) is fulfilled, the cooling process 
acts too quickly after the initial collision to allow the main 
shock going far enough from the white dwarf surface to fulfil the 
first criterion (see Eq.(p^). Also, it should be noted that when 
both cyclotron and bremsstrahlung are effective, considering the 
stabilizing effect of cyclotron losses, the secondary shock does 
not appear for a high magnetic field for the same reason. 

From the simulations, the detailed X-ray spectrum emitted 
by the PSAC over an oscillation can be computed and is shown 
in Fig. The spectrum evolves strongly through the different 
phases of the oscillation. When only the primary shock is 
present {t = 6.35 s), the emitted spectrum is softer than when 
the secondary shock has appeared {t = 6A- 6.43 s). The relative 
difference between the mean spectrum calculated over a couple 
of oscillations and the spectrum emitted by the theoretical 
stationary column is presented at the bottom of Fig. The 
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Fig. 4. Evolution with time of the position of the main shock (blue filled 
points), the secondary shock (red filled points) and the rarefaction wave 
(grey dashed line) in the same condition as in Fig.[^ The Mach number 
value for the two shocks are also shown by the open symbols. Different 
critical times through the oscillation are labelled by letters and discussed 
in the text. 


difference between the two spectra is quite negligible at low 
energy, but at high energy, the mean spectrum is significantly 
lower than the stationary one. Thus, if we attempt to evaluate 
the white dwarf mass through the mean spectrum, we expect to 
find a lower mass. 



E,[eV] 


Fig. 5. Top: Evolution of the emitted bremsstrahlung spectrum through 
the oscillation for the same parameters and the times indicated in the 
inset of Eig.|^ The mean spectrum calculated over a couple of oscilla¬ 
tions is represented by the thick grey line. Bottom: Relative difference 
between the mean spectrum and the spectrum emitted by the theoretical 
stationary column. 


Considering both cooling instabilities, the column structure, 
when cooled mainly by the bremsstrahlung process, is quite dif¬ 
ferent from the steady state solution. This strongly affects both 
light curves and emitted spectra. 


4. Influence of system parameters on luminosity 
variations 

In Paper I, observational X-ray light curves of a set of po¬ 
lars were analysed to search for oscillations as predicted by 
the cooling instability model. To directly compare these obser¬ 
vational constraints with the theoretical expectations from the 
shock model, we used the Hades simulations to compute the am¬ 
plitudes and frequencies of the expected oscillations more pre¬ 
cisely for different sets of parameters that are relevant to the po¬ 
lars, considering both bremsstrahlung and cyclotron cooling. 

4.1. Bremsstrahlung dominated columns 

The temporal characteristics of the front shock position and 
the computed PSAC X-ray flux emerging from the simulations 
were analysed using fast Fourier transforms (FFT) to better 
understand the properties of the quasi-periodic phenomena. 
Representative power spectra of the shock front position and 
X-ray luminosity are presented in Fig. where we consider a 
white dwarf of 0.8 M© with an accretion rate of m = 10 g/cm^s 
and only taking the bremsstrahlung losses into account. The 
FFT covers simulation data obtained at 0.01 s resolution over 
2048 time steps at long times between t = 19.76 s and r = 30 s. 

Several discrete frequencies emerge from theses analyses, 
and the same frequencies are present in both signals but with 
different relative amplitudes. The different amplitudes are 
explained by considering that not only does the variation in the 
emission volume modify the luminosity but also the physical 
conditions of temperature and density throughout the post-shock 
region. For instance, the post-shock temperature depends on the 
pre-shock velocity in the shock frame. Thus, when the shock 
is going upstream (resp. downstream), the infalling flow goes 
faster (resp. slower) than the free-fall velocity, and the post¬ 
shock temperature is higher (resp. lower). This phenomenon is 
clearly visible in the temperature profiles presented in Fig. 
Amongst the most significant frequencies, the lowest is found 
at 2.24 Hz (period of 0.45 s) (mode/). This frequency is very 
visible in the first instants (see Fig.|^ and remains present later. 
It is correlated with the setting up of the PSAC structure: when 
the flow impacts the surface, the shock goes away from the 
surface and progressively slows down. At ^ = 0.24 s, the shock 
falls back and reaches the minimum position at r = 0.42 s before 
bouncing against the secondary shock. The second strongest 
mode is present at a frequency of 6.7 Hz (period of 0.15 s) (mode 
a). This typical timescale corresponds to the cooling timescale 
that is calculated using Eq. 0- Modes of higher frequencies 
are also present and can make a significant contribution to the 
oscillations. A mode (b) has been emphasized in Fig. since 
it has a strong contribution for higher magnetic field strength. 
For the regime with pure bremsstrahlung cooling, the FFT 
exhibits the same overall behaviour for the different white dwarf 
masses tested (Mwd = 0.4,0.6,0.8 M©). However, the relative 
contribution of each mode can be slightly different, and it can be 
demonstrated how the frequency depends on the accretion rate 
and the white dwarf mass. 
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Fig. 6. Fast Fourier transform of the front shock position (left) and the 
X-ray luminosity (right) respectively normalized by their mean values 
for the same parameters as in Fig. The frequencies are identical for 
both signals but with different amplitudes. The main modes discussed 
in the text are referred as (f), (a), and (b). 
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Considering the bremsstrahlung cooling as dominant (see 
Regime 1, Fig.[^, the mechanism that induces the oscillations 
is the cooling instability with the oscillation frequency expected 
to evolve as the inverse of the cooling time (see Eq. (|^) 

(18) 


Introducing the expression of the free-fall velocity coupled to 
the white dwarf mass-radius relation ( |Nauenber^ 1 197 2| ) , the 
frequency depends on the specific accretion rate and the white 
dwarf mass as 


y = A^fn 


1.44MoV^^ 

Mwd / 


1.44Mo 

Mwd 


4 / 3 ' 


ll/2 


(19) 


where Aq is a constant. In Fig.[^ we represent the dominant os¬ 
cillation frequency of the X-ray luminosity (mode a), extracted 
from simulations for three different white dwarf masses and 
accretion rates. The dependence from Eq. ( p^ is plotted as solid 
lines and gives a good estimate of the oscillation frequency. The 
constant Aq, which provides the best fit, is Aq = 0.9,0.9, and 
0.7 cmVg for m = 0.1,1, and 10 g/cm^/s, respectively, while the 
expected theoretical value is Ao,r/z = 0.5 cm^g. These values 
are slightly different since the theoretical value is evaluated just 
behind the shock and does not take the strong variations along 
the column into account. 


These results imply that, when QPOs are observed, their 
frequencies should allow constraining the value of the accretion 
rate and white dwarf mass (see Paper I). 


4.2. PSAC efficiently cooled by the cyclotron process 

When the magnetic field is strong enough, the cyclotron 
cooling becomes efficient and a two-process cooling has to be 
considered. The losses for the bremsstrahlung and cyclotron 
processes are presented in Pig. as a function of the position 
in the column for a typical configuration. The cyclotron losses 
are higher than the bremsstrahlung ones at the shock front in 
this case since 6s = 1.2, which corresponds to a magnetic field 
of 30.3 MG, assuming a column section of 10^^ cm^. Because 
the bremsstrahlung mechanism is more sensitive to variation in 
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Fig. 7. Variation in the main oscillation frequency (mode a) of the X-ray 
light curve as a function of the white dwarf mass, for different accretion 
rates in the case of pure bremsstrahlung cooling. The curves represent 
the best fits following the theoretical tendency given by the relation ( p^ . 


density than temperature (see Eq. ([^), the losses progressively 
increase from the shock to the surface of the white dwarf along 
with the density. In contrast, since cyclotron mechanism is more 
sensitive to variations in temperature than in density (see Eq. 
([^), it is more efficient at the shock where matter is hotter than 
closer to the white dwarf surface. Thus, 6s, defined as the ratio 
of cooling times (see Eq. ([^), is always maximum at the shock 
and then decreases along the column towards the white dwarf 
surface. 



Normalised position in the column 

Fig. 8. Variation in the bremsstrahlung (solid blue line) and cyclotron 
(dotted red line) losses along the column assuming a white dwarf of 
0.8 Mq, an accretion rate of m = 10 g/cm^s and 6s = 1.2 (corresponding 
to a magnetic field of 30.3 MG) at t = 100 s. The abscisses have been 
normalized to the front shock position (Xs ~ 7.7 X 10^ cm). 


The analysis in Paper I gives upper limits for the amplitude 
of quasi-periodic oscillations in X-rays. Up to now, oscillations 
in light curves have only been observed in the optical range for 
a few selected sources. In the context of the shock instability 
model for polars, it should be emphasized that optical variations 
associated to cyclotron are expected together with the X-ray 
ones linked to bremsstrahlung. From the simulations, we 
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extracted both components separately and compare them here in 
terms of amplitude and frequency. 

In Fig. the luminosity in the optical range due to the 
cyclotron losses is compared to the X-ray emission due to 
the bremsstrahlung process, considering a white dwarf with 
Mwd = 0-8 Mq and an accretion rate of 10 g/cm^s, for a 
typical magnetic field of 18.7 MG (6s = 0.3). The behaviour 
in the two energy ranges is very similar. We present in Fig. [T^ 
the corresponding Fourier transforms of the cyclotron and 
bremsstrahlung emissions, where the power shown is the 
oscillation amplitude normalized by the mean luminosity. They 
present the exact same frequencies of oscillations but with an 
oscillation amplitude in the dominant mode (a), which is lower 
for the cyclotron signal (~ 18%) than the bremsstrahlung one 
(~ 28%). This behaviour has been systematically observed in 
all simulations. Thus, we expect to observe stronger oscillations 
in the X- ray light curve th an in the optical domain. As already 
noted by |Wu et al. | ( |1992| ), we confirm that the two signals are 
slightly out of synchronization. 



time (s) 


Fig. 9. Variation in the (0.5-10 keV) bremsstrahlung (upper blue line) 
and cyclotron (lower red line) luminosity considering a white dwarf of 
0.8 Mq, an accretion rate of m = 10 g/cm^/s, and a magnetic field of 
18.7 MG (6s = 0.3). 
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Fig. 10. Fourier transform of the (0.5-10 keV) bremsstrahlung luminos¬ 
ity and cyclotron total luminosity with the same parameters as in Fig.[^ 


Considering the perturbative analysis, the system is expected 
to stabilize as cyclotron becomes a strong contribution in the 
radiati ve losses ( [Saxton et al. |1998t|Saxton |1999||Saxton & Wu| 
[]|1999| ). In the results presented in this work, for a white dwarf 
of 0.8 Mo, the oscillations have been considered as completely 
damped when the front shock oscillation is only a few meshes 
in length. It corresponds to a luminosity variation of less than 
~ 0.1%. We observed that for the accretion rate of m ~ 1 
g/cm^s, the oscillations are effectively damped for 6s ~ 1.5; and 
for m ~ 10 g/cm^/s, the oscillations are almost damped for the 
same value of 6s. 

As the magnetic field strength increases, the low-order 
modes are expected to be rapidly damped. This type of be¬ 
haviour is found in our numerical simulations. The typical 
evolution of the X-ray light curve FFT is presented in Fig.[TT]for 
various values of 6s that correspond to a magnetic field strength 
ranging from 0 to 28 MG, considering an Mwd = 0.8 M© 
white dwarf and an accretion rate of m = 10 g/cm^/s. For weak 
magnetic fields, the dominant mode is around (y ^ 6-1 Hz), but 
for higher fields, this mode decreases strongly and disappears, 
while higher modes (v ^ 20 - 25 Hz) appear. 



Frequency (Hz) 

Fig. 11. Changes in the FFT of the X-ray light curve for a growing 
magnetic field characterized by 6s, for a white dwarf of 0.8 M© and an 
accretion rate of m = 10 g/cm^s. The range in 6s from 0 to 1 corre¬ 
sponds to a magnetic field from 0 to 28 MG. The different frequency 
modes are identified as (a), (b), and (c) and listed in Table|^ 


In Fig. [m it is clear that the distribution of the various 
modes and their identification are not obvious. However, a few 
main modes can be tracked when the magnetic field grows. The 
frequencies and relative amplitudes of the three main modes 
are reported in Tables and and plotted in Fig. [T^ We see 
that Mode (a) decreases monotonically as the magnetic field 
increases. Its amplitude becomes negligible for a magnetic field 
around 23 MG. Mode (b) increases initially and then drops 
suddenly around the same critical magnetic field as Mode (a). 
Mode (c) has the same behaviour as Mode (b), but its amplitude 
keeps increasing when the first two modes decrease strongly. 
Its amplitude ultimately decreases for higher magnetic field 
amplitude. The different behaviours exhibited by these first 
three modes is interesting since it reveals the great richness 
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of the non-linear dynamics of the PS AC. Moreover, it should 
be emphasized that progressively, as the magnetic field grows, 
the power is distributed in higher order modes. It should be 
noted that this non-linear behaviour is coherent with the linear 
stability studies published by |Saxton et al. ( 1998| ). In Tables 
and the total relative amplitude is also given as a quadratic 
sum of the amplitudes in the different modes when the power is 
split into different frequencies. 



Fig. 12. Variation in the relative amplitude of the X-ray oscillations for 
the three main modes observed in the FFT for a white dwarf of 0.8 M© 
and an accretion rate of m = 10 g/cm^/s. The total relative amplitude is 
given by the full black line. 


For a given mode, a frequency variation with the magnetic 
field strength is observed in the simulations that can be explained 
analytically. Under the assumption that oscillations are induced 
by the cooling instability, the expected frequency is directly con¬ 
nected to the cooling time by the relation 

y oc l/4ooi = -U + U = [1 + f(5)] ^ (20) 

%rem ^cyc %rem 

where 6(B) is the ratio between the bremsstrahlung and cyclotron 
cooling time in each spatial position. If we suppose that the typ¬ 
ical timescales calculated behind the shock are representative of 
the value along the accretion column, then 6(B) ^ 6s, and the 
relation can be written as 


^ ~ ^brem,s 


1 -h 


( B 
110 MG/ 


( 21 ) 


where ybrem,s is the oscillation frequency at the front shock in the 
asymptotic bremsstrahlung dominated regime, and Q is a con¬ 
stant that depends on the accretion rate and the free fall velocity 
as 


Co ^ 2 X 10 


,-4 


r 1 

-17/40 

fn 

_ 10^^ cm^. 


.1 g/cmVs. 


Vff 


10^ cm/s 


-37/20 


117/20 


( 22 ) 


The specific behaviour of the oscillation mode fr equencies 
predicted here is different to the results presented b y |Wu et al. 
n (119921 ) in the context of a stochastic accretion. In |Wu et al.~ 
( |l992| ), the suggested simplified evolution law v cx (1 +CofiO is 


Article number, page 10 of 13 



Fig. 13. Evolution of oscillation frequency with magnetic field for the 
first two modes (a and b) for a white dwarf of 0.8 M© and an accretion 
rate of m = 10 g/cm^/s. The lines are the best fits according to the 
predicted dependency of Eq. ( |^ . 


only based on phenomenological assumptions and not justified 
by physical considerations. 


The evolution of the frequency for a range of magnetic fields 
from 0 to 22 MG is presented in Fig. 13 for the two main modes 
(a) and (b), assuming a 0.8 M© white dwarf and an accretion 
rate of 10 g/cm^/s. It should be noted that for strongest magnetic 
fields, the amplitude of those modes decreases drastically so no 
oscillation frequency can be found. The curves in Fig. [^repre¬ 
sent the best fits using the analytical trend presented in Eq. ^T\) . 
The constant determined from these fits is Co = 0.022 and is the 
same for the two modes. In this specific configuration, the ex¬ 
pected analytical value of the constant as derived from Eq. ( j^ 
is Co,th = 0.061. The difference with the fitted value from the nu¬ 
merical simulations comes from the structure of the post-shock 
region. Indeed, to theoretically infer the characteristic timescale, 
only the values just behind the shock have been taken into ac¬ 
count (see also 4.1). 


5. Conclusions 

In this work, a numerical study of the dynamics of PSAC has 
been presented. Radiation hydrodynamic simulations have 
shown that the coupling between magnetic fields, radiation 
processes, and extreme hydrodynamics leads to a complex 
dynamics and a wide variety of dynamical behaviours. To 
validate the numerical simulations, they have been compared to 
theoretical stability criteria, analytical results, and dimensional 
tendencies. Good agreement has been obtained, in general. To 
link simulation results to observational data presented in Paper I, 
a post-treatment of numerical results were performed to extract 
and study the X-ray and optical luminosities and their dynamics, 
the X-ray spectral luminosity, and the frequency spectrum. 
Given the observational data (see Paper I), we have focused 
our analysis on the two most important regimes of PSAC: the 
bremsstrahlung-dominated regime (weak magnetic field and/or 
high accretion rate) and the coupled bremsstrahlung-cyclotron 
regime (strong magnetic fields and/or lower accretion rate). 

In the bremsstrahlung-dominated regime, one of the most 
important results is the development of the catastrophic in¬ 
stability, which leads to the formation of secondary shocks. 
According to radiation hydrodynamic simulations, those addi- 
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tional shocks modify the dynamics of the primary shock, as well 
as density and temperature profiles. It induces modifications 
on the total and spectral emissivities and plays an important 
role in the oscillations of the shock height. The complex cou¬ 
pling between the secondary shock dynamics and the cooling 
instability development leads to quasi-periodic oscillations in 
the light curves. Specifically, the secondary shock implies a 
harder X-ray spectrum than in the stationary regime, which can 
affect the mass determination from the classical X-ray approach. 
The presence and role of the secondary shocks in PSAC have 
not been studied much as yet. In principle, the presence of 
the secondary shock can also be revealed through variations 
induced in the X-ray spectrum across the oscillations. However, 
the sensitivity of the present X-ray telescopes is not yet high 
enough to detect such variations. Laboratory experiments, using 
powerful lasers already in use or near completion such as NIF 
or LMJ, are a promising way to directly investigate these shocks. 


regimes, the plasma parameter (noted p = PthIPmag, where Pth 
and Pmag are respectively the thermal and magnetic pressure) is 
significa ntly different. It is in th e range y^xTauri ~ 10“^ - 10^ for 
T Tauri (Matsakos et al. ||2013), but only ySpoiar ~ 10“"^ - 10“^ 
for polars. therefore the magnetic collimation is much more 
efficient for polars, and 2D effects are expected to be less 
important. In T Tauri, low p values have been shown to lead 
to a fragmentation of the column into independent oscillating 
tubes or fibrils ( [Orlando et al. ||2010| [Matsakos et al. ||2013| ). If 
out of phase, the oscillations of the different fibrils can cancel 
each other out. Such an effect can also occur in polar PSAC and 
would also induce a frequency widening, as well as a significant 
decrease in the overall QPO amplitudes. These effects will be 
the subject of an upcoming study where numerical simulations 
will take multi-dimensional effects into account, as well as 
non-homogeneous accretion and a more detailed treatment of 
the boundary conditions at the white dwarf photosphere. 


In the coupled bremsstrahlung-cyclotron regime, numerical 
simulations show the stabilizing effect of the cyclotron process. 
The predicted amplitude of oscillations in the X-ray light curve 
considering a pure bremsstrahlung regime is ~ 30%. As the 
magnetic field grows, the amplitude of oscillation decreases and 
a critical magnetic field depending on the accretion parameters 
can be found where no oscillation persists. Thanks to numerical 
simulations, a more detailed study of the distribution of power 
in the different modes of oscillations has been explored. In 
particular, the low-order modes are the first to be damped as 
the magnetic field grows, and the power is then distributed in 
higher order modes. These two effects have been theoretically 
expected in the linear regime and persist with non-linear effects. 
An important correlation between the magnetic field and the 
oscillation frequency has been established. A comparison 
between the numerical results and the theoretical relations based 
on dimensional arguments leads to very good agreement. 

Numerical simulations exhibit strong oscillations in the 
X-ray luminosity with the same frequencies as in the optical 
luminosity curve. They show the presence of QPOs with a wider 
range of frequencies in the X-ray domain (0.6-50 Hz) than com¬ 
monly assumed. This wider range, which has never been fully 
explored observationally until now, has been investigated in 
Paper 1 with negative results. There is therefore no compelling 
evidence yet of significant X-ray oscillations, in contrast to what 
is expected from shock instabilities. These numerical results 
raise important questions on the origin of physical processes 
that lead to QPO. 


The inconsistency between X-ray and optical observational 
data, on the one hand, and between numerical expectations and 
observations, on the other, highlights that the physics of QPO 
is still not completely understood. The situation is somewhat 
similar among classical T Tauri stars for which shock-heated 


accreting material is expected to give rise to QPOs ( [Koldoba 
et al. [[200^ [Sacco et al. ||2008) that are not yet observed 


( [Drake et al. [|20Q9t [Gunther et [|20TQ| . In this context, 2D 
axisymmetric MHD simulations have shown that significant 
distortions of the accretion structure can occur as a result of 
the gas pressure. At the bottom of the column, accumulation 
of matter may lead to outflows parallel to the star surface that 
can increase the cooling, and the oscillations are expected to 
be reduced or suppressed (Orlando 2010, 2013). However, the 
analogy between polars and T Tauri stars is somewhat limited. 
Indeed, for the typical magnetic field and velocity and density 


Additionally, it should be emphasized that the PSAC 
presents interesting similarity properties that allow us to repro¬ 
duce a scaled model in laboratory on measured scales ( [Falize et 


al. |201 lb] [Busschaert et al. |2013b[ ). This is the main ob jective 


of a laboratory astrophysics project: the POLAR project (Falize 
et al. |20lT^ . The possibility of producing relevant experiments 


has b een demonstrated on the LULI2000 facility ( [Falize et~S7 


20121. The experimental way is a complementary approach to 


increase our understanding of the radiation hydrodynamics of 
accretion shock. Thus the QPO physics will be explored using 
dedicated designed experiments with powerful lasers. 
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Table 1. Characteristics of the main oscillation modes extracted from the simulations of an accretion column onto a 0.8 M© white dwarf with an 
accretion rate m - I g/cm^s and for different magnetic field strengths. The first column gives the value of the ratio 6s defined in iD. the second 
column gives the corresponding magnetic field strength considering an accretion surface of S = 10^^ cm^, and the third column gives the value 
of the theoretical stationary shock height. The subsequent columns give for the first three main modes (a, b, and c), the frequency and the relative 
amplitude of the oscillations in the (0.5-10 keV) X-ray range (Lx) and the optical cyclotron flux (Lc). The last two columns give the quadratic 
mean value of the relative amplitude integrated over all the oscillation modes, including those not reported in the table. 
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Table 2. Same as Table [^except for m = 10 g/cm^/s 
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